Topological aspects of quantum spin Hall effect in graphene: 
Z 2 topological order and spin Chern number 
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For generic time-reversal invariant systems with spin-orbit couplings, we clarify a close relationship 
between the Z2 topological order and the spin Chern number proposed by Kane and Mele and 
by Sheng et at, respectively, in the quantum spin Hall effect. It turns out that a global gauge 
transformation connects different spin Chern numbers (even integers) modulo 4, which implies that 
the spin Chern number and the Z2 topological order yield the same classification. We present a 
method of computing spin Chern numbers and demonstrate it in single and double plane of graphene. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.61.Wp, 85.75.-d 



Topological orders 0, play a crucial role in the clas- 
sification of various phases in low dimensional systems. 
The integer quantum Hall effect (IQHE) is one of the 
most typical examples 0, , in which the quantized Hall 
conductance is given by a topological invariant, Chern 
number, due to the Berry phase [5j induced in the Bril- 
louin zone [fj . Such a topological feature should be more 
fundamental, since it has a close relationship with the 
parity anomaly of Dirac fermions [zj 0, 0. 

Recently, the spin Hall effect EjjEl El El has been 
attracting much current interest as a new device of so- 
called spintronics. In particular, Kane and Mele 0,0] 
have found a new class of insulator showing the quantum 
spin Hall (QSH) effect ES E3, E3 

which should be re- 
alized in graphene with spin-orbit couplings. They have 
pointed out |l5j that the QSH state can be specified by a 
Z2 topological order which is inherent in time-reversal 
(T) invariant systems. This study is of fundamental 
importance, since the Z2 order is involved with the Z2 
anomaly of Major ana fermions |l9U2fl|. 

On the other hand, Sheng et al. |21j have recently 
computed the spin Hall conductance by imposing spin- 
dependent twisted boundary condition, generalizing the 
idea of Niu et al. 0)113] • They have shown that it is given 
by a Chern number which is referred to as a spin Chern 
number below. This seems very natural, since the QSH 
effect is a spin-related version of IQHE. The spin Chern 
number for graphene computed by Sheng et al. indeed 
has a good correspondence with the classification by Z 2 . 

However, the Chern number is specified by the set of 
integers Z, not by Z 2 . Although the studies by Sheng et 
al. [13] suggest a close relationship between two topo- 
logical orders, natural questions arise: How does the con- 
cept Z2 enter into the classification by Chern numbers, or 
otherwise, does the spin Chern number carry additional 
information? 

In this letter, we clarify the relationship for generic T 
invariant systems. We show that while two sectors in 



the Z2 classification are separated by topological changes 
due to bulk gap-closing phenomena, each of these sectors 
is further divided into many sectors by boundary-induced 
topological changes in the spin Chern number classifica- 
tion. The latter is an artifact which is due to calcula- 
tions in finite size systems with broken translational in- 
variance. Therefore, the different spin Chern numbers in 
each sector of Z2 describe the same topologically ordered 
states of the bulk. 

Consider generic electron systems on a lattice with T 
symmetry, described by the Hamiltonian H . Denote the 
electron creation operator at jth site as cj = (cjjjcj^). 
Then, T transformation is defined by Cj — ► Tcj with 
T = ia 2 K, where the Pauli matrix a 2 operates the spin 
space and K stands for the complex conjugation opera- 
tor. Let H(k) be the Fourier-transformed Hamiltonian 
defined by H = £ fc J(k)H(k)c(k) and let |n(fc)) be 
an eigenstate of H(k). Assume that the ground state 
is composed of an M-dimensional multiplet of degen- 
erate single particle states which is a generalized non- 
interacting Fermi sea Q ■ Kane and Mele E3 have found 
that the T invariant systems have two kinds of impor- 
tant states belonging to "even" subspace and "odd" sub- 
space: The states in the even subspace have the prop- 
erty that \n(k)) and T\n(k)} are identical, which occurs 
when TTi(k)T~ 1 = TL(k). By definition, the states at 
k = always belong to this subspace. The odd subspace 
has the property that the multiplet \n(k)) are orthogo- 
nal to the multiplet T\n(k)) . These special subspaces can 
be detected by the pfaffian pkm(&) = pf (n(k)\T\m(k)) . 
Namely, pkm(^) = 1 in the even subspace and in the 
odd subspace. Kane and Mele have claimed that the 
number of zeros of PKu(k) which always appear as T 
pairs ±fc* with opposite vorticities is a topological invari- 
ant for T invariant systems. Specifically, if the number 
of zeros in half the Brillouin zone is 1 (0) mod 2, the 
ground state is in the QSH (insulating) phase. 

We now turn to the spin Chern number proposed by 
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Sheng et al. |2fJ. According to their formulation, we 
impose spin-dependent (-independent) twisted boundary 
condition along l(2)-direction 
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(1) 



where a set of integers j = (j'1,.7'2) specifies the site 
and f and 2 stand for the unit vectors in 1- and 2- 
directions, respectively. Let Tt(9) denote the twisted 
Hamiltonian and let \n{9)) be corresponding eigenstate. 
It follows from Eq. (JJJ that T transformation induces 
TU{e u e 2 )T- 1 = H(6i,-9 2 ) and therefore we can al- 
ways choose \n{0 1,-62)) = T\n(9i,9 2 )) except for 62 = 
0. The states on the line 92 — belong to the even 
subspace. Below, the torus spanned by 9 is referred to 
as twist space. The boundary condition JQ) enables us 
to define the spin Chern number, but the cost we have 
to pay is the broken translational invariance. In other 
words, we slightly break T invariance at the boundary, 
which leads to nontrivial spin Chern numbers. The av- 
erage over the twist angles recovers T invariance. 

Let us define a pfaffian for the present twisted system 
as a function of the twist angles: 

p(6) = pi(n(6)\T\m(0)), n,m = 1,- ■ ■ ,M, (2) 

where M is a number of one particle states below the 
Fermi energy with the gap opening condition. Since the 
line 92 = belongs to the even subspace, the zeros of the 
pfaffian in the 62 > or < twist space can move only 
among the same half twist space keeping their vorticities, 
and never cross the 62 — line. Therefore, one T pair of 
zeros are never annihilated, like those of the KM pfaffian 
PkmQ^). Furthermore, the two zeros in the same half 
twist space are never annihilated unless they have the 
opposite vorticity. This is a crucial difference between 
the KM pfaffian and the twist pfaffian. The number of 
zeros in the twist pfaffian should be classified by even 
integers (half of them, by Z) . 

Are these pfaffians topologically different quantities? 
The answer is no. To show this, let us consider a global 
gauge transformation Cj — > g(ip)cj, where 



9{<f) 



ia ip 



cos^ + itr sm<p. 



(3) 



This transformation replaces the the Pauli matrices in 
the spin-orbit couplings into g t (<p)crg(tp) = (cos 2ip a 1 — 
sin2(/?er 3 , a 2 , cos2ipa 3 + sin 2<,5 c 1 ). On the other hand, 
since Eq. © is an orthogonal transformation, one- 
parameter family of transformed Hamiltonian, denoted 
by H v or W{6) below, is equivalent. The pfaffian © is 
also invariant. Therefore, when we are interested in the 
bulk properties, we can deal with any Hamiltonian H v . 

So far we have discussed the bulk properties. However, 
if we consider finite periodic systems like Eq. QJ , a family 
of the Hamiltonian H v behaves as different models. It 
follows from Eq. |Q that the gauge transformation © 



is commutative with e 2 , but not with e l6ia . This tells 
us that spin-dependent twisted boundary condition is not 
invariant under the gauge transformation @ , and breaks 
the gauge-equivalence of the Hamiltonian H v which the 
bulk systems should have. 

To understand this, the following alternative consider- 
ation may be useful: If we want to study bulk properties 
of T invariant systems, we can start with any of H v . For 
one H v with ip fixed, let us impose the twisted boundary 
condition After that, we can make a gauge trans- 
formation © back to H°. Then, we can deal with the 
original Hamiltonian H , but with a gauge-dependent 
twisted boundary condition for x-direction; 



iOi (cos 2tp(T 3 — sin 2ipa 1 ) 



(4) 



Namely, the gauge equivalence is broken only by the 
boundary condition in 1-direction. Now, imagine a situa- 
tion that at ip = the pfaffian (J2J) has one T pair of zeros. 
We denote them as (6l,±8 2 ) with vorticity ±m. Let us 
change ip smoothly from to ir/2. Then, it follows from 
Eq. (0J that at ip — tt/2 the coordinate of the torus is 
changed from (6\, 6*2) into (—6*1, #2) and therefore, we find 
the zeros at (— 9\, ±6*2") with vorticity =Fto. We thus have 
a mapping (9{,±9*) -> (-f?*,^) for tp = -> tt/2. 
Namely, the zero in the #2 > (< 0) space moves into 
the 02 < (> 0) space, and thus the zeros can move 
in the whole twist space like those of the KM pfaffian. 
During the mapping, there should occur a topological 
change, but it is attributed to the boundary, i.e, an ar- 
tifact of broken translational invariance, and half of the 
zeros of twist pfaffian should be also classified by Z 2 , by 
taking into account the gauge transformation 

Using this pfaffian, we next show that its zeros can be 
detectable by computing the spin Chern number which 
seems much easier to perform than searching them di- 
rectly. The spin Chern number \2l\ is defined by c s = 
^-J d 2 9F 12 {9), where F 12 (9) = d l A 2 {6) - 8 2 A 1 {6) is 
the field strength due to the U(l) part of the (non- 
Abelian) Berry potential Q, A M (8) = tr(n(6>)|<9 M |m(6»)) 
with dfj, = d/89^. First, we will show the relation- 
ship between the zeros of the twist pfaffian and the spin 
Chern number c s . For the time being, we fix p. The de- 
generate ground state as the M-dimensional multiplct 
has a local U(M) gauge degree of freedom, \n{6)) — > 
Em \ m (9))Vmn{9) , where V(9) is a unitary matrix 0. 
Let us denote V(9) = e la ^/ M V{9), where dct V(fc) = f. 
This transformation induces A^(9) — > A fJ/ (6)+ 18^(6) to 
the Berry potential. If one can make gauge- fixing globally 
over the whole twist space, the Chern number is proved 
to be zero. Only if the global gauge-fixing is impossible, 
the Chern number can be nonzero. 

Among various kinds of gauge-fixing, we can use 
the gauge that p{9) is real positive, because p{9) — > 
p(6) det V{6) = p(8)e la ^\ This rule can fix the gauge 
of the Berry potential except for p(6) — 0. Therefore, 
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nontrivial spin Chern number is due to an obstruction to 
the smooth gauge-fixing by the twist pfaffian. This cor- 
respondence also proves that the spin Chern number is 
an even integer, since the pfaffian @ always has T pair 
zeros, and since Fi2(0i,— 62) = ^12(61, 82)- 

Let us now change p. At ip = 0, we obtain a certain 
integral c s . Remember that at <p = 7r/2, the coordinate of 
the torus is changed into (—6*1, 6*2). Therefore, we have a 
mapping c s — > — c s for (61,62) — > (—61,62)- As in the case 
of the pfaffian , we expect topological changes during 
the mapping. However, as stressed, these changes are ac- 
companied by no gap-closing in the bulk spectrum, just 
induced by the symmetry breaking boundary term which 
is an artifact to define the Chern numbers. Therefore, 
the states with ±c s should belong to an equivalent topo- 
logical sector. Since the minimum nonzero spin Chern 
number is 2, we expect c s mod 4 (if we define the spin 
Chern number in half the twist space, mod 2) classifies 
the topological sectors. 

So far we have discussed the Z2 characteristics of the 
spin Chern number c s . We next present several exam- 
ples. To this end, we employ an efficient method of com- 
puting Chern numbers proposed in Ref. [23j. We first 
discretize the twist space [0, 2ir] ® [0, 2tt] into a square 
lattice such that M = 2-Kj tl /N^, where — 1, ••■ , 
[24]]. We denote the sites on this lattice as 81 with 
£ = 1,2, ••• ,N 1 N 2 . We next define a U(l) link vari- 
able associated with the ground state multiplet of the 
dimension M, 

U^8 e ) = \detU fi (6 i )\- 1 detU^e), 

where L^(^) m „ = (m(6i)\n(6 e + ft)) with n,m = 
1, • • • ,M denotes the (non-Abelian) Berry link variable. 
Here, fx is the vector in li direction with — 2ir/N fl . 
Next define the lattice field strength, 

F 12 (6 e ) = \nUx(6 t )U2(6 t + l)U^(8 e + 2)U^(6 e ), 

where we choose the branch of the logarithm as 
1-^12(0^)1 < 71 '■ Finally, manifestly gauge invariant lat- 
tice Chern number is obtained: 

Cs = ^E F ^). (5) 

t 

As shown in [23j |. the spin Chern number thus defined 
is strictly integral. To see this, let us introduce a lattice 
gauge potential A ll (8g) = In Up (6 e) which is also defined 
in |Ap(0^)| < 7T. Note that this field is periodic, A fl (di + 
Np) = A fl (8 e ). Then, we readily find 

F 12 (6 t ) = AiA 2 (0/) - A 2 A 1 (6 ( ) + 2win 12 (8t), (6) 

where A p stands for the difference operator and ni 2 (6^) 
is a local integral field which is referred to as n-field. 
Finally, we reach c s = ^^12(81). This completes the 
proof that the spin Chern number is integral. While the 



n-field depends on a gauge we adopt, the sum is invariant. 
For the T invariant systems, the pfaffian ^) is very useful 
for the gauge-fixing also for the lattice computation. In 
the continuum theory, we have stressed that the zeros of 
the pfaffian play a central role in the Z2 classification. 
Since such zeros occur at several specific points in the 
twist space, it is very hard to search them numerically. 

Contrary to this feature in the continuum approach, 
we can detect the zeros in the lattice approach as fol- 
lows: Suppose that we obtain an exact Chern number 
using Eq. (J5J with sufficiently large 7V M . Since the lat- 
tice Chern number is topological (integral), which implies 
that even if we slightly change the lattice (e.g, change the 
lattice size or infinitesimally shift the lattice), the Chern 
number remains unchanged. Next, let us compute the 
n-field in the gauge that the pfaffian is real positive. If 
the zeros of the pfaffian happen to locate on sites of the 
present lattice, we cannot make gauge-fixing. Even in 
such cases, we can always avoid the zeros of the pfaf- 
fian by redefining the lattice with the spin Chern num- 
ber kept unchanged. Thus we can always compute the 
well-defined n-field configuration. If the Chern number 
is nonzero, there exist nonzero n-field somewhere in the 
twist space. These nonzero n-field occur in general near 
the zeros of the pfaffian: Thus, without searching zeros 
of pfaffian in the continuum twist space, we can observe 
the trace of such zeros as nonzero n-fields. 

Now let us study a graphene model with spin-orbit 
couplings Hi HQ; 

h = - 1 4°3 + ~h v *° Yl c i CT ■ ( dk i x dik ^ c i 

(id) V «M» 

+ iVr ^ 4 (a- x dij ) 3 Cj +v s ^2 s g n J C J C J ' ( 7 ) 

(i,j) 3 

where cj = ( c jj? c jj) is the electron creation opera- 
tor at site j on the honeycomb lattice, sgnj denotes 
1 (—1) if j belongs to the A (B) sublattice, and e£y 
stands for the vector from j to i sites. The first term 
is the nearest neighbor hopping, the second term is the 
s 2 -conserving next-nearest neighbor spin-orbit coupling, 
whereas the third term is the Rashba spin-orbit cou- 
pling, and final term is a on-site staggered pot ential. 
Analyzing the KM pfaffian, Kane and Mele jl4] have 
derived the phase diagram: It is in the QSH phase for 
16^3^0 — v s \ > yjvl + 9 Vp and in the insulating phase 
otherwise. Sheng et. al. [2l| have computed the spin 
Chern number c s = 2 in the QSH phase and in the 
insulating phase. 

For numerical computations, it is convenient to use 
the momentum &2 0, 113 instead of 02 because of the 
translational invariance along this direction even with the 
boundary condition First, we show the spectrum in 
Fig. [H at 0i = as a function of k 2 - The left belongs 
to the QSH phase with c s = 2, whereas the right to the 
insulating phase with c s = 0. This topological change is 
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FIG. 1: Spectrum for 6\ = as a function of k-z. Parameters 
used are V so = O.lt, v s = 0.3t, and Vr = O.lt (left), Vr = 
0.225207t (middle), and Vr = 0.3t (right). 



due to the gap-closing in the bulk spectrum, as shown in 
the middle in Fig. ^ Therefore, the phase with c s = 2 is 
topologically distinguishable from the phase with c s — 0. 
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FIG. 2: Left: Spectrum for nonzero <p = 7r/4 as a function 
of k-i at Q\ — it/2. Other parameters are the same as those 
of the left in Fig. Right two: the n-field configuration 
corresponding to the left in Fig. Q The left (right) is at tp = 
(tt/2). The white (black) circle denotes n = 1 ( — 1), while 
the blank means n = 0. We used the meshes N\ = N2 = 10. 



Contrary to this, how about the phases with c s = ±2? 
As we have mentioned, the phase with c s changes into — c s 
when we vary ip from to tt/2. In this process, boundary- 
induced topological changes must occur. We show in Fig. 
121 the spectrum cut at 6\ —tt/2 for 99 = 7t/4. We indeed 
observe a gap-closing at finite 9\, and the spin Chern 
number c s = 2 for < <p < 7r/4 is changed into c s = — 2 
for 7r/4 < (p < tt/2. As stressed, this change is attributed 
to the boundary (edge states in Fig. [2J, not to the bulk, 
and we conclude that the phase c s = ±2 is classified as 
the same QSH phase. These spin Chern numbers c s = ±2 
are well visible by the n-field. In Fig. [2J we also show 
the n-field for tp = and 7r/2 cases in the QSH phase. 
The points of nonzero n-field is closely related with the 
positions of the pfaffian zeros. We also note that in Eq. 
ijSJl net contributions to the nonzero Chern number is just 
from two points. 

Next, let us study a bilayer graphene. Suppose that 
we have two decoupled sheets of graphene described by 
H Vi with i = 1,2 whose lattices include A, B sites and 
A, B sites, respectively. For simplicity, we take into ac- 
count only the interlayer coupling 71 between A and B 
0: Via = 7i J2j c |l j C2B <i + k- c -> wnere » = 1, 2 in c l:j 
indicate the ith sheet. Now make the gauge transforma- 
tion © separately for each sheet to obtain the same H° 



as Eq. (JJJ. Then, we have identical bilayer system H° ® 
H° coupled by V 12 = 71 J2 3 ; c l^ j 9( l Pi)9 t ( l P2)c2Bd+h..c. 
with two independent boundary conditions c ij+ii j = 
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FIG. 3: The n-field configuration for 71 = O.lf. Other pa- 
rameters are same as those of the left in Fig. Q Left: (pi = 
and y?2 = 0. Middle: ipi = n/2 and <p2 = 0. Right: ipi — -k/4 
and ip2 = — tt/4. In the case ipi = n/2 and <f2 = tt/2 we have 
the same figure as the left but with black circles. We used the 
meshes Ni = N 2 = 10. 

For the same parameters as those of the left in Fig. ^ 
the spin Chern number is, of course, Cg = 2 + 2 = 4 in 
the limit 71 = 0. This spin Chern number remains un- 
changed for small but finite interlayer coupling 71 . How- 
ever, taking into account the gauge transformation g(tfi), 
the spin Chern number changes. In Fig. |3| we show ex- 
amples of the n-field for 71 = O.lt. We have the spin 
Chern numbers 2 + 2 = 4, 2-2 = 0, and -2-2 = -4: All 
of them are denoted as c s = mod 4, which belong to the 
insulating phase. Detail analysis of this model including 
the interlayer coupling 73 will be published elsewhere. 

Finally, we comment that the QSH effect is understood 
by the edge states and therefore, it is interesting to 
establish the bulk-edge correspondence for T invari- 
ant systems with respect to Z2 . We also mention that Fu 
and Kane j2^ and Moore and Balents [2t| have recently 
discussed the relationship between the Z2 order and the 
spin Chern number, and reached a similar conclusion. 
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